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I. INTRODUCTION 

This paper deals with the calculation of the spectral zeta function 

oo 

n=l ™ 

associated with the eigenvalues E n of the inhomogeneous Hclmholtz equation on a d-dimensional domain SI 

(—A d )ip n (xi,...,Xd) = E n S{xi,...,x < i)ip n (xi,...,Xd) , (2) 

where = jp- + ' • "\"m~ is the laplacian operator in d-dimensions and ip n (xi, ■ ■ ■ , %d) are its eigenfunctions (the index 

n refers to the entire set of quantum numbers that specify a solution.). E(zi, . . . , Xd) is a density (E(a;i, . . . , x d ) > 
for x G 

In general, spectral zeta functions can be built using the eigenvalues associated to a given Schrodinger operator: 
for example Voros has studied the zeta function of anharmonic oscillators 0, [|| , Steiner has obtained sum rules for 
confinement potentials of the form V(r) — gr p , with p > and g > [H, Berry has studied the zeta function of 
Aharonov-Bohm quantum billiards at integer values of s extending the results of Itzykson, Moussa and LuckQ, 
who had derived an explicit expression for the sum of inverse powers of the Laplacian on a simply connected region 
with Dirichlet boundary conditions at the border; this problem has been later studied by Steiner [6J, obtaining exact 
expressions for the spectral zeta function at positive integer values, and by Elizalde, Leseduarte and Romeo[7[, who 
derived exact sum rules for the zeros of the Bessel functions J„ using a recursive formula (these sum rules were already 
obtained by in Q, considering an infinite square well of radius R); Steiner has also obtained exact expressions for the 
Selberg's zeta function for compact Riemann surfac es p (in particular the case of the two-torus is discussed in detail 
in the review article by Arendt et al.); Crandall[Iol| has calculated the exact values of the spectral zeta functions 
at integer values for several one dimensional potentials H; Kvitsinsky [l2| has obtained sum rules for two-dimensional 
domains which are close to the unit disc; Mezincescu (l3j ]. and later Bender and Wang[T3|. have used the spectral zeta 
functions to obtain strong evidence of the reality of the spectrum of certain 'PT-symmetric hamiltonians. 

Spectral zeta functions are also useful tools in the calculation of the Casimir effect, since the Casimir energy is 
obtained as 1/2 Z(— 1/2). The series given in eq.([T]) diverges in this case and one needs to analytically continue it 
to negative values of s to obtain a finite physical value for the Casimir energy: this is the essence of zeta function 
regularization method [l5l - fl9j . Applications of these techniques may be found in [20i - [23| . 

Unfortunately, the direct use of eq.([T|) is possible only if the eigenvalues E n are known exactly; in this paper we 
develop a perturbative approach to the calculation of Z(s) that allows one to obtain explicitly the corrections to the 
spectral zeta function to a given order in the perturbation, bypassing eq. |T]). In our calculation we assume that Jl is 
a d-cube with a density E which is a perturbation of a uniform density (in two dimensions may also be obtained 
from the small deformation of a square, a circle or a domain where a basis is available). 

We have derived a general expression for the spectral zeta function for problems with small inhomogeneities and/or 
deformations, which we have worked out explicitly to second order in the perturbation. We have applied our formula 
to different problems, obtaining in many cases the analytic continuation to all values of s (this has been generally 
accomplished only to first order, and to second order in a specific example a one dimensional string). For the case of 
a string with piecewise constant density, we have reproduced to first order in the perturbation the results obtained 
by Hadasz, Lambiase and Nesterenko[2J|. Although most calculations of Casimir energies for inhomogeneous systems 
concern the special case of piecewise constant densities [25l - |28j . our formalism applies to more general densities. 
In particular, we have observed that the Casimir energy of an inhomogeneous string, calculated to first order in 
perturbation theory, in some cases contains irremovable divergences, when either Dirichlet or Neumann boundary 
conditions are used, but that these singularities cancel in the sum of the two zeta functions corresponding to the 
different boundary conditions. 

This paper is organized as follows: in section [II] we briefly introduce the Rayleigh-Schrodinger perturbation theory 
for a general hermitian operator; in section Hill we obtain an explicit formula for the spectral zeta function using the 
perturbative expansion of section [TTJ in section HVl we discuss several applications in one, two and three dimensions; 
finally, in section [V] we draw our conclusions. 



In particular he has showed that Z(l) for the power-law oscillator obtained by Voros holds for arbitrary power-laws with positive real 
exponents. See also ref. [Tl| . 
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II. RAYLEIGH-SCHRODINGER PERTURBATION THEORY 

Raylcigh-Schrodinger perturbation theory is a standard tool of quantum mechanics: it allows to obtain perturbative 
expressions for the eigenvalues and the eigenstates of a hermitian operator C = A+rjB, defined as the sum of a solvable 
hermitian operator A and a of perturbation rjB (\rj\ <c 1), which typically does not commute with A. 

In this framework, using the basis of the eigenstates of A, one obtains an explicit expression for the eigenvalues of 
C, as a power series in 77, which to second order reads 



, ... . 2 v-^ \{n\B\k)\ 2 



k^n 



The formula above assumes that \n) is a non-degenerate eigenstate of A. 

Assuming that for an arbitrary small (but finite) 77, rjB is a perturbation for all the states belonging to the spectrum 
of C, we may write the trace of C s as a power series in 77 using the explicit perturbative expression for c n : 



Tr C s = 



< + V sa^{n\B\n) + ^^—^-a'- 2 (n\B\n) a +*f £E ^ \{n\B\k)? + O W] 

n k^n 



(4) 



where the last term has been symmetrized with respect to the summation indices. Notice that this term is well defined 
even when |fc) and |n) are degenerate states, even though the corresponding expression in eq.Q diverges in this case. 
In a similar fashion one may extend these results to the operator f(C) = f(A + r/B), obtaining 



Tr /(C) = £ 



(5) 



f(a n ) + V f'(a n ){n\B\n) + ^/"( a „)(n|i3|n) 2 + ^ £ £ ~ M _ 7 ~ + ■ ■ 

n k^tn 

III. SPECTRAL ZETA FUNCTIONS 

We work with the spectral zeta function of Eq.([T]) associated with the eigenvalues E n of the hermitian operator 
O = ~^^(~ on ^ ne d-dimensional cube ft (Ad is the d-dimensional Laplacian and S(a^i, . . . ,Xd) is a positive 

function on Q). The spectrum of O is bounded from below, with strictly positive eigenvalues for Dirichlet boundary 
conditions for (xi,X2, ■ ■ ■ , Xd) G <9f2. 

If £ is a physical density, then the eigensolutions of O provide the normal modes of an inhomogeneous d-dimensional 
cube; in two dimensions, S may also be obtained by conformally mapping the eigenvalue problem defined on an 
arbitrary two dimensional domain onto the square. 

The specific form of O given above is obtained from the Helmholtz equation ([2]) after recasting it into the equivalent 
form 

1 1 

—=(-A d )—=4> n (xi,. ..,x d )= E n <p n (xi, ...,x d ) , (6) 
v S v S 

where cf) n (xi, . . . , Xd) = VS ip n (x\, . . . , Xd). The advantage of working with eq.([B]) rather than with eq. resides in 
the manifestly hermitian form of the operator O = ^^( — ^rf)"^- 
Alternatively, one may obtain a still equivalent equation 

1 1 1 

= S(xi, . . ■ ,X d ) r = £ n {Xi, ... ,X d ) = —£, n {xi,...,X d ) , (7) 



where £„(xi, . . . , x d ) = y/ (-Ad) ip n (xi, . . .,x d ). 

In this case one works with the hermitian operator Q = . 1 S(xi, . . . , x d ) , 1 El- Notice that the eigenvalues 

of eq. (J7|) are the reciprocal of the eigenvalues of eqs.Q and (|6]). 



Another possibility would be working with the operator O 1 = vE(-Aj) 1 v / S, which being the inverse of O is isospectral to Q. 
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In d dimensions, the asymptotic behavior of the energies of eq. ([2]) for n ^> 1 is described by Weyl's law, E n w 
47r(r(l + d/2) f) 2/rf , where n is the number of states with energy smaller than E n . As a result, the series X^^Li "W 
will converge to a finite value only for s > d/2. 

Using the invariance of the trace with respect to unitary transformations, for s > d/2, we may express the spectral 
zeta function directly in terms of the diagonal matrix elements of the operators 0~ s and Q s : 

Z(s)= 2_j {ni,...,n d \6~ s \ni,...,n d ) = E (ni, . . . ,n d \Q s \ni, . . . ,n d ) , s > d/2 (8) 
m,...,nd ni,...,na 

where \n±, . . . , n d ) is an eigenstate of (—A d ) on the d-cube ft and e ni ,...,n d the corresponding eigenvalue (when possible 
we will use the compact notation |n) and e„ to refer to \n\, . . . ,n d ) and to c ni nd ). For integer values of s, with 
s > d/2, eq.© provides exact sum rules. 

We will now discuss the case in which S(xi, . . . , x d ) = 1 + 5Yi(x\, . . . , x d ), with \5T,\ <C 1, corresponding to a small 
perturbation of a uniform problem; taking into account this fact we write: 



{-A d ) ^=35 V 3 ^ 

and treat the second term in the parenthesis as a perturbation. 

After identifying A = (-Ad)" 1 and B = (-A d ,)- 1 / 2 S'E(-A d )- 1 / 2 we then obtain 

Z(s) = TrQ s = £ ( S - f E E 4 ~': C "~ a HgS|fc) 2 + O (9) 

where O [<^ 3 ] refers only to terms containing off-diagonal matrix elements, since the purely diagonal contributions 
have been included to all orders. We define the first term in the expression above 

Z (diag)^ / MEM 



n 

to be the component of Z(s) which only contains diagonal matrix elements. 

There is an important aspect concerning eq.©: since the matrix elements of (n\B\k) are always finite in the limit 
k, n — > oo 0, the perturbative hierarchy in the expression for Z(s) is not destroyed by an unbounded growth of the 
perturbative corrections with respect to the quantum numbers and therefore eq.([9|) provides a genuine perturbative 
expansion in this case. 

On the other hand, we can easily find an example where the perturbative expansion for the spectral zeta function 

-2 , -2 . 

breaks down: in the case of the dimensional anharmonic oscillator, with hamiltonian H — p ^ — h Arc (0 < A <C 1), 
the perturbative corrections to the eigenvalues of H are polynomials in the quantum number n of increasing order 

x o 3 3\ , 2 / 17 , 51 2 59 21 \ 
E n = (n + l/2) + X^-n 2 + -n+-j + A 2 {-—n 3 - -n 2 - —n - — J + . . . (10) 

and therefore for any arbitrary small A the expansion breaks down for sufficiently large n. Therefore, for any finite 
value of A one cannot apply perturbation theory to all the spectrum of H . Notice however that our formalism applies 

-2 

to the hamiltonian H — + r\ V(x) (0 < rj <C 1), where 



V(x) = 



v{x^j X. m i n ^ X ^ %max 
OO X <C X m i n , X > &rn.ax 



and v(x) is a potential bounded from below. In this case the asymptotic behavior of the spectrum is determined by 
the confinement of the particle between the infinite walls and does not depend on the perturbation. This property is 
sufficient to make rjV(x) a perturbation for all the spectrum, for a given infinitesimal rj. 
Let us now consider the operator: 

1 



W^Q- l = ^K^j- rA/FAd), (11) 

2-i\X\, . . . ,X d ) 



It is easy to check that lim n ->oo("|ffi|n) = ^O d ^{ x )d d x and limn—xxj (n|<5£|fc) = for a fixed k, with k ^ n; in ref. l29t explicit 

expressions for these matrix elements for Dirichlet be in one dimension are derived. 
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which is isospectral to — "^( — ^d)^ used in ref.[3C|. The heat kernel of W is defined as 

#(*) = £« 



-X„t 



(12) 



where A„ are the eigenvalues of W. 
In order to apply eq.(|5]) we write 



W = (-A d ) + v/Fa^j 



1 



and identify A = (-A d ) and _B = yJ(-A d ) 
Using eq.© the heat kernel reads 



-(e n +(n\B\n)) 



E(xi,...,x d ) 
t 



. . . ,Xd) 



v/FaF 



n fc^n 



{n\B\kf 



Notice that for t — > + we have 



EE 



(13) 



using the completeness of the basis; in this limit we have the perturbative expansion for the heat kernel: 



A-(i)*5>" 



l-t(n|S|n) + —(n\B 2 \n) 



The spectral zeta functions associated with the eigenvalues of W may be obtained from the heat kernel as 

poo 



Z( S ) 



n») Jo 



t*- 1 K(t) dt 



(14) 



(15) 



(16) 



(17) 



IV. APPLICATIONS 

In this section we consider several applications of the formulas obtained in the previous sections. 



A. One dimension: Inhomogeneous strings 

Before discussing specific examples of inhomogeneous strings, it is useful to make some general considerations: the 
Casimir energy of an inhomogeneous string is obtained calculating the value of its spectral zeta function at s = —1/2. 
Assuming that the density of the string is a perturbation of a constant density (£(#) »1 + <5E) and working to first 
order in perturbation theory, we have Z(s) s» 1+s ^"J as ^ , where e„ are the eigenvalues of the homogeneous string 

subject to specific boundary conditions and \n) are its eigenstates. 

The analytic continuation of this series allows one to obtain the Casimir energy of the string, directly evaluating 
Z{— 1/2): however, if Z[s) contains terms which are proportional to £(2s+2), the corresponding Z{— 1/2), calculated to 
first order, diverges, making the physical interpretation problematic. Previous examples of Casimir energies containing 
irremovable divergences are known: for example, Sen calculated the Casimir energy of a circle in two dimensions, 
finding that its Casimir energy is infinite [3ll [32| ; Bender and Milton have calculated the Casimir energy of a massless 
scalar field in a hyperspherical shell in D spatial dimensions and observed that it diverges when D was a positive even 
integer [33j. 

It is easy to see that for an inhomogeneous string, to first order in the perturbative expansion a divergent Casimir 
energy is obtained when the matrix element (n|<5£|n) contains a term going as 1/n 2 . In the case of Dirichlet boundary 
conditions, this matrix element may be expressed in terms of the series [291 ] 



(n!5S|n) 



2L 



8Y,(x)dx 



fe=0 



(27m) 



(18) 



6 



where 5Y,( 2k+lS> (x) = d d tk+T ; therefore the Casimir energy of an inhomogeneous string is finite only if 8T![x)\ x —l — 
8Y,' (x)\ x= -l. This condition is fulfilled if the perturbation is either an odd function of x or if <5E(x) is flat at the 
border: 8Y,'{x)\ x = L = 6Z>{x)\ x =- L = 0. 

Therefore the Casimir energy of the inhomogeneous string contains an irremovable divergence when none of these 
conditions is met; on the other hand, irremovable divergences may still show up to higher orders, even when these are 
absent to first order, as we will see in an example. 



1. Casimir energy of a string with piecewise constant density 



The Casimir energy of a string with piecewise constant density and four different sets of boundary conditions at 
its ends (Dirichlet-Dirichlet, Neumann-Neumann, Dirichlet-Neumann and Neumann-Dirichlct) has been calculated in 
ref . [24| ; the density of the string is 



E(aO = I Pl ' ~ L < X ~ X T ° 
{ P2 , x < x < L 



To uniform our notation to that of ref.[2J|, we define pi } 2 = l/uf 2> R = r = xq + R/2 and E = oj 2 . 
a, (3 and Sv are defined as in ref. [24j : 

Vi — V2 r R— r r R — r 
ov — , a — 1 , p — . 

V\ + V2 V\ V2 V\ V2 

For Dirichlet boundary conditions, the exact eigenvalues are the solutions of the trascendental equation 

sin au + Sv sin /3a; = , 

which corresponds to the first of eqns. (2.4) of ref. [23 |. 

For a string of arbitrary density, the energies go asymptotically as E„ » ' wnere &{L) = f_ L \J ' p{x)dx (see 
for instance [34{): in the present case <j{L) — a and therefore E n ps n 2 7r 2 /a 2 for n — > oo, which is the behavior of the 
spectrum of a uniform string with constant density S = CT ^ = 

We may therefore write 

S(x) w t + SY,(x) 

and apply our formulas to first order in <5£, for which only diagonal matrix elements of the perturbation are needed. 
For Dirichlet boundary conditions, the basis is 

(D) , , IT . nir(x + L) 

w„ '(x) = \ — sm , 

Vn y ' V L 2L 

with n = 1,2,.... 

Correspondingly, the eigenvalues of the negative laplacian with Dirichlet be are 

2 2 
(DD) _ n n 

4L 2 ' 

The matrix elements of S in this basis are 

TTn(vl(L-r)+vl(L+r))+L(v 1 -v 2 )(vi+V2)sm{ 7 ' n{L L + r] ) _ , 



lEln'' 



,(DD) 



( < ,?-t,I)((nW)»in( -^+-»"-" , > )+(n'-n)sm( ^ +r »" + "'> )) 



n 7^ n 1 



Notice that (nlEIn)^ 00 -* does not contain terms going like 1/n 2 for n — > oo and therefore we expect that the Casimir 
energy of this string will be finite. 

A simple calculation, carried out to order 5v, yields 



,(DD) 



ir((2s) +i s Sv (Li 2s+ i (e £ ^ \ - Li 2s+ i (e ( ^ Vj 



7r a 
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where e„ = n 2 7r 2 /4L 2 and Li t/ (z) = X^fcli i s tne polylogarithmic function. 

We observe that this equation provides an analytic continuation of the series to negative values of s; in particular, 
for s = —1/2 this formula returns the Casimir energy of the string, which reads 



i _ Svttm(^) 

yfDD) I i /r,\ _ i V 2a 7 



E (DD) = V(DD) ( _ 1/2) = ___ (ig) 

The same result can be obtained regularizing the series, evaluated at s = —1/2, with a cutoff function, f(a,n) 
e~ an , and taking the limit a — > + at the end of the calculation: 



E, 



(DD) _ 

2a 2 a 4a 24a 



l-i(n| [E/E-l] |n)| (DD) 



a->0 H 



Notice that the finite part of this expression agrees with the result of eq.(|T9|. 

We can compare this result with the analogous result of eq.(3.5) and (3.6) of ref.[24]; taking Sv <C 1 and calculating 
(3.5) to order Sv we obtain the result of our eq. tjlU)) . 

We may also calculate the exact value of the spectral zeta function for this string at s = 1: 

W , DD V 1N I (a 2 - /3 2 Sv 2 ) ,. (a + (3Sv) 2 / / \ / i(2+aw*\\ 

z (1) ° g iiF - 6d - + ,fo J. - 1)» ( L - ( e "~ )- L,3 ( e ^ ))' 

This expression holds to all orders in <5f. 

The calculation for the cases of Neumann-Neumann, Neumann-Dirichlet, Dirichlet- Neumann and periodic-periodic 
boundary conditions can be done in a similar fashion. Here we just report the results. 

The spectral zeta functions for these cases, calculated to first order, are 

Z^( S ) = TT- 2s a 2s (^s) + i5v{2K)- 2s - l S e-^a 2s ($ [e~^ ,2s + 1, 

_ e <• $ e = , 2s + 1, - I - e °> Li 2s+ i « J + e <« Li 2s+ i « j I , 
Z< DN '(s) = (4 s - l)7r- 2s a 2;5 C(2s) 

o ^ of HLA ( if>JL l\ Mt ( iSjL 1 

+ (5-dtt 2s_1 sa 2s ( e _ ~$ ( -e «, 2s + 1, - 1 + e~$ I -e~ , 2s + 1, - 

Z< pp )( S ) = 2 1 - 2s vr- 2s a 2s C(2 S ) , 

where $ (z, s, a) = 53fcLo ( a +fc) s is the Lerch transcendent function. The case of Neumann-Dirichlet be is simply 
obtained from the case Dirichlet-Neumann be with the substitution Sv — > —Sv. 
The Casimir energies for these cases are then 



4 NN) 



P (DN) 



E (PP) 



tt ^an(ff) 

24a 4a 

48a 4a 

faBec(jg) 

48a ' 4a 

7T 

6a 



These results agree to order Sv with the expressions calculated in ref. [24] , with the exception of the PP case which 
was not studied there. 
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Exact sum rules can also be obtained using eq. 

„2 _ «2r,,2\ f„ , ax.,\2 



z (nn) x _ v ; - ifo-i ' ' Li 3 (e stptj — - Li 3 (e — — 



(<fc 2 - 1) 

7 (dn) n \ _ (a - 8Sv)(a + 8Sv) iSv(a + fiSv) 2 

U 2 - 2^2 + ^3 (^2 _ !)2 



z( pp )(i) 

Notice that 



i7T(aV+l)( Q +g) / i(Q+g)(fe+l)^ 1 

e 2 (°+- 3 *-) $ e °+' 3a « .3, - 
V 2 



a 2 - /3 2 (5t/ 
12 - 12<Su 2 



. u + + / i(a + g)(i^ + l)ir 1 



Z (PP)(1) = I( Z (DD) (1) + Z (NN) (1) ^ 



3. yl slightly inhomogeneous string 
We consider a string of length 2L (\x\ < L) and with density 

7TX 

= 1 + r]6T, = 1 + rysin — , 

where |ry| < 1 (this condition is needed to enforce S(x) > for x £ (—L,L)). 

In this case the matrix elements of ST, in the basis of the homogeneous string are particularly simple 

(n\5T\m) = — - [5 mt n+i + 8 m , n -i] , 

and allow one to calculate exactly the value of the spectral zeta function of the string at the first few positive integers 
values of s: 

Z(l) = \l 2 

, , 8L 4 8n 2 L 4 24n 2 L A 
Z{2) = 1 ' 



Z{3) 



45 3tt 2 tt 4 
64L 6 16 V 2 L 6 64r] 2 L 6 720t? 2 L 6 



945 15tt 2 tt 4 tt 6 

Using our eq.([5]) we may calculate the spectral zeta function of this string to second order: 

, w _ , (£) - m + a)*± + o M ^ 



2n + 1 

n— 1 



Let us define: 



* (i H«)-E^-f:^c(, + i-a)4 + f:^ (ccj+i-»)-d 

n=l j=0 j'=0 



and 



1 00 1 

*(a) = *«(a) - *( 2 )(a) = ^ ( C (2j + 2 - a) - 1) . (20) 



3 2 2 J+! 

Notice that ^(a) is singular at odd values of a, a = 1, 3, 5, . . . and that ^(2) = \&(4) = • • • = 1. 
Thus the spectral zeta function reads 

j \ 2s / T \ 2s 



Z{s)=A s - C(2s) +^s4 s - i - #(2 -2s). 
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In particular, around s = — 1/2 we have: 

z <») ~ ^ Trrm - m + Hi ( 3 * ( t) + * - 31 ) + ° (<s + 1/2)2) ■ 

where 7 ss 0.577216 is the Euler-Mascheroni constant. 

In this case a divergence in the Casimir energy shows up to second order in perturbation theory, since STf {x)\ x —l = 
ST/{x)\ x ——i,. Clearly the divergences which emerge at a given order in the perturbative expansion cannot be canceled 
by divergent contributions stemming from higher orders, given the different dependence on the expansion parameter. 
The cancellation of these divergences therefore requires considering equally divergent terms which may show up in 
the spectral zeta function corresponding to different boundary conditions evaluated to the same order in perturbation 
theory. 

We will now prove that the sum Z( DD )(s) + Z< NN )(s) does not contain a singularity at s = —1/2 to first order in 
perturbation theory @. To prove this point consider 

z {DD+NN) {s) = Z (DD) (s) + Z (NN) (s) 

1 (n\ [E/E- 1] |n)( DD ) 



8- 



" { \ ' " ) I e " DD) ) 



2 



(n,u\ [E/E-l] 



\(NN) 



.(NN)V' 
~-n,u J 



/ (NN)V' 
1 tn.u I 



where 



p (nn) = J t?~ , " = 0,1,2,... , u=l 
n ' u \ , n = l,2,... , u = 2 

are the eigenvalues of the negative ID laplacian with Neumann-Neumann boundary conditions and E = er(L) 2 /(4L 2 ). 
Notice that 

(DD) _ (NN) (DD) _ (NN) 

t 2n — e n,l > fc 2n-l — fc n,2 



Moreover 



where 



L ' 



1 



n = , u = 1 
^("OH ^ cos2 F - n = l,2,... , u = l , 

1 • (In — 1)ttx in o 

Tf sm 2L > n = l,2,... , m = 2 

are the eigenfunctions of the negative ID laplacian with Neumann-Neumann boundary conditions. 

Using these identities we may write the spectral zeta function for the one-dimensional string to first order as 

z iDD +NN){s) w 2t s (l + JL( L [E(x)/S - 1] dx) £ * , (21) 

where X^^Li (ddi- — vr _2s (2L) 2s C(2s) is the spectral zeta function of the homogeneous string. Therefore the only 
singularity of the expression in eq. (|21[) is the singularity of the Riemann zeta function at s = 1/2 and the corresponding 
Casimir energy is always finite, as anticipated. 



4 In two and three dimensions this particular combination of boundary conditions is needed to evaluate the Casimir energy of the 
electromagnetic field. 



10 



*(a) 



3.0 




1 






2.5 










2.0 


- 


\ 






1.5 
1.0 
0.5 


': 
















-2 

-0.5 




2 


4 


6 


-1.0 












FIG 


. 1. *(a) of eg. (|20l). 



a 



3. Borg string 



We consider a string of unit length (L = 1/2) with density 



E(x) 



(l + a(x+l)) 



where \x\ < 1/2 and Dirichlet boundary conditions at its ends, a is a free parameter (a > —1). It has been proved 
by Borg that this string is isospectral to a homogeneous string of the same length and unit density [35j : therefore the 
exact energies are 



E n = 7r 2 n 2 



and the corresponding spectral zeta functions are easily calculated: 

C(2s) 



Z 



Borg 



00 = 



r 2.s 



It is possible to obtain an explicit expression for the matrix elements of E in the basis with Dirichlet boundary 
conditions; for instance, the diagonal matrix elements read 



(nlSln) = 



27r 2 (q + l)n 2 
3a 4 



a z + 27r(a + l)n Ci 



where Ci(a;) and Si(x) are the cosine and sine integrals. 
Using these matrix elements we obtain the identity 



cos 



2(q + l)nn 



/2nn 
V a 



Ci 



2nn 



sm 



^ OO 

Z B org{l) = ~ = 



(n|E|n> 



6 

OO 

E 



2(a + 1) 
3a 4 



o + 27r(a + l)n Ci 



2(a + l)ri7r 



-Ci 



2n7r 
a 



27rn 



si(^)-siP (a + 1)n7r 



27rn 
a 



which we have numerically verified with high accuracy. 
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The isospectrality implies that the spectral zeta function of the Borg string is independent of a: using this property 
and working to a given order in a, we may obtain non trivial mathematical relations. For example, if we select the 
contribution of order a 2 in the perturbative expression for the spectral zeta function we obtain the identity: 

~ k 2 n 2 ((-l) fe +" - 1) (n 2 - 2s - k 2 - 2s ) 7T 4 , N 5tt 2 , , 

E(s) = V ^ J -^-= '- = C(2s) (2s + l , (22) 

W ^ 2(k 2 -n 2 ) 5 768^ ' 256^ 1 " ' V ' 

which provides the analytic continuation of the double series to negative values of s (see Figj2]) . 

We may also obtain non trivial identities working with the heat kernel: in this case we need to evaluate the matrix 



^m-nts..l n-=V(-;E*)5VHE*): 

3q 4 (a 2 +3a+3)a 2 tt 2 (a 4 +5a 3 + 10a 2 + 10a+5)n 2 _ 

. _ . 2^(a + l)V ~ (a+1) 2 H 5(a+l) 2 ' " " ™ 

n\W\m) - ^ 2 56« 2 TO 4 n 4 (l2a 2 (ro 2 4V)+(«+l)(7r 2 («+l) 2 ^ > 

If we use eq. (fH)) . to order a 2 we obtain the identity: 

"'" (( -" ~ e } - 4 *• I* (»■ -° 



n— 2 m— 1 



(24) 



where $3(2, g) = X)^L-oo 9™ e 2mz is the Jacobi theta function. We have tested this identity numerically for several 
values of t. 

Using eq. flTTl) inside eq. (j2"4")l we have 

y y ^ — — [ 4 — - — ui((2 8 )-i((2( >+ i)), (25) 

h k ^=i 3(k 2 -k 2 ) 5 768 U ' 256 U { " { ' 

After combining this identity with eg. (1221) we obtain an alternative series representation for the Riemann zeta 
function: 

00 n— 1 



c( S )^ee uj ; — — > . (26) 

n n=2k=l (k 2 -n 2 ) 
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4- A string with oscillating density 



We consider now a string with rapidly oscillating density 

= 2 + t) sin (27r(a; + L)/e) , 

where |x| < L and e — » + . 

This particular problem has been studied by Castro and Zuazua in ref . [36[ and more recently by myself in ref. [2!| . 
The solutions to this problem with wavelength com par able to the typical size of the density oscillations are localized at 
the ends of the string |3(i - l38j . In particular, in ref. |29l | we have obtained both numerical and analytic approximations 
to the low part of the spectrum of this string (the analytic result improves the results previously obtained by Castro 
and Zuazua in ref. [36|). 

In this case the non-zero matrix elements of the density are 

2- ir?(-l) n sin(^|pl , n = m,n=l/e 

i ^ e 3 rin 2 sinf ) sinf ) , w , 

'nl^Un. I 2 + — 1 -—^J—±iJ- , n = m,n^l/e 

8e 3 ?)TO7i(-l)"' + "(cos( 7r( ^+^ ) ) — cos(j (m+ra)) cos( ""^r )) . , s , .„<_ 

where e = e/2L. 

As in the example of the piecewise constant string we may calculate 



2Le , - 

<t(L) = + 2 

7T 



, ((s + 2)L-2£)tt 2 v \ E ((eL-2(L + £))TT f 2r, 



AsL V + 2 J \ 4eL 1 rj + 2 



where E (</>|m) is the incomplete elliptic integral of second kind. 

The asymptotic behavior of the spectrum of this string is E n w n 2 ir 2 /a(L) 2 , which is also the spectrum of a 
homogeneous string with density £ = a{L) 2 j AL 2 . 

Working to first order in rj we have 

n=l { £ " £ ™ J 

*E((^J + ^+1 Sm liJ Sln (,ilJ 

n— 1 v \ / 

This expression contains the series 

$ (^>EirT. ( 27 ) 

which we may express as 



/ J e^-n 2 — 1 



00 2— 2s 

rr * 
e 2 n 2 — 1 

n=l 



We] 2-2s 



* — ' e z rr, z — 1 ^ — ' £^n z — 1 

71=1 n=[l/e] + l 



where [a] is the integer part of a. 

(2) 

The function $s (e) may be cast in the form 



* (2 W) = E EW = Ei^cf2 (fc + S ) 



n=[l/e]+l *;=0 fe=0 



where C(^ a ) is the Hurwitz zeta function. Notice that $( 2 )(s, e) has an infinite number of singularities, located at 
semi- integer values s = i — j, with j = 0, 1, . . . 
For | ry | <C 1 the spectral zeta of this string is 

8L 2 V _ , ^ 2 3s-i L 2s (lxl 



„, . XL, -\ >( . . £ u ^° "Zv" . /7T\ . / 7T£ \ , 
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FIG. 3. $(s, e) °f eq. (f27)l as function of s for e = 2/101; the horizontal line is lim s _nx> $ ( '{s,e) 



For s — > —1/2, this spectral zeta function to first order in rj reads 



Z{s) 



2AV2L 



£ 3 r]sm (|) sin (g^) 
8\/2L 



oo , 

$W(-l/2, e -)+ £ 3S=aC(2*-l, 





1" 




^2fc- 1, 











W 2s + 1 ' 



where H n = X)fe=i \ is the harmonic number of order n. Notice that the term corresponding to k = 1 in the series is 
treated separately because of the singularity of the Hurwitz zeta function. 
It is also useful to write: 



$«(-l/2, £ -)-ff w =£ 



n 2 e — 1 n 



Hn] i H_i Hi 



2e 2 



2e 2 



2e 2 2e 2 



Since Z{— 1/2) diverges for n ^ 1/e or for £/eX non integer, the Casimir energy of this string is infinite in these 
cases. Moreover, in this example, there is an infinite sequence of singularities of the spectral zeta function located at 
negative semi-integer values of s. 

It is easy to see that this divergence is absent if the original density is replaced by a piecewise constant curve, for 
instance £(x) = T,( x ^ h+1 ) for Xk < x < Xk+i- In this case the diagonal matrix elements of do not contain 
the 1/n 2 term and the resulting spectral functions arc finite at s = —1/2. Calling h = x^+i — X). = 2L/N , this result 
holds for arbitrary small (but finite) h (or equivalently for arbitrary large N). Here h acts as a cutoff of the highly 
excited states with quantum numbers n 2§> l/h. 

Let us discuss in detail this approach; in this case we have 



a(L) = 2V2L + ^L sin (|) cse sin (g) + O ( V 



and £ = a(L) 2 /4L 2 . 

The diagonal matrix elements of £ may be evaluated explicitly and read: 



(n|E|n) 



T ? sin(|)sin(M) ( 27mC sc (^) - TV sin (^) (esc (2=fa) + esc (^f±^))) 



2itnN 
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Working to first order in 77 we have 

*(•>»>-'*-*■ (?)«(£) 



sin I 5 ) ) C(2s) 



-2s-l 



71=1 



. /7rn 
sin ^— l l esc 



7r — 7ren\ / 7ren + 7r 



The function W(n) = sin (^) (csc (^^p) +csc(^±zi)) i s periodic with period TV, W(n + AT) = W(ra), and 
therefore we may write 



Z W « 2- (2 + £ S sin (J) esc sin (g) ) C (2a) 



+ 77 (-2 3s - 2 ) tt" 25 -^ sin ( ^ ) L 2s sin ( ^ iV 



*r-2a-l 



AT-1 



"Hat) l csc (^v-J 



esc 



For a finite N, this expression is finite at s = —1/2: the divergent behavior of the Casimir energy of the string with 
a smooth density is recovered only in the limit N — > 00, where the sum over n diverges. 

As a further example, we may calculate the spectral sum rule at s = 1 for a string with an arbitrary density of 
period A and length 2L = 1: 



E(*) = ?+E 



27T7X , . 27T7X 

a, cos — t 1- b n sin — - — 



where 



M v-^ a 



a,jA . 2Ljn 



U jL7T 

and M = J^ L Y,(x)dx is the total mass of the string. 

The diagonal matrix elements of the density in this case are 



sin ■ 



00 A 3 n 2 a ? sin ( 
1 ivi \ a ° V A / 

(n|E|n) = o--2^ 



and 



n— 1 ' 



2 ^87rj 3 i 3 -27rA 2 jin 2 
A 3 a, sin A 2 a 3 cos (^) ALaj sin 



47T 3 j 3 L 



27T 2 j 2 



37TJ 



Notice that Z(l) reduces to the sum rule of a homogeneous string for A = 2L/N with N integer. 

B. Inhomogeneous drums 

In this section we will discuss the calculation of the spectral zeta functions of inhomogeneous drums. 

1. Spectral zeta function of a square membrane 

We consider a square of side 2L: in this case the eigenvalues of the negative laplacian (with Dirichlet boundary 
conditions at the border) are known exactly and read 



t 1 2 \ 
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with n x ,n y — 1,2,.... 

The series defining the spectral zeta function for this problem, 



oo 1 

*°M= E is- 



E s 

n x ,n H = l ""> n s 

may be analytically continued to all s. Ziff (39j | and Steiner Q report the spectral zeta function of the square: using 
the notation of Ziff it reads 

Zn{s)= ^) 2S (/3( S )C( S )-C(2 S )) , 

where (3(s) = X)^Lo( — l) n (2n + l) _s - Notice that s = 1/2 and s = 1 are the only singularities of Zq(s). 
Around s = 1 the zeta function behaves as 

.An f 1 

where An = AL 2 is the area of the square. 7 ss 0.577216 is the Euler-Mascheroni gamma and 71(a) is the Stieltjes 
constant defined as 71(a) = lim s _>i (jry — C( s i a )) (C( s i a ) — J2^=o ( n +a) s * s ^ ne Hurwitz zeta function). 



2. Spectral zeta function on an arbitrary simply connected domain in two dimensions 



Let T> be an arbitrary simply connected domain of the plane and call E n the eigenvalues of the Helmholtz equation on 
this domain, assuming Dirichlet boundary conditions on dV. Riemann's theorem grants the existence of a conformal 
map from T> to any other simply connected region of the plane and in particular to {x, y) € f2, O = {|x| < L, \y\ < L}: 

w = u + iv = f(z) , 

(u,v) e V. 

Under this map, the original Helmholtz equation on T> transforms into the Helmholtz equation on f2, in presence 

2 

of a variable density y) = 4^ (see ref. 30]). We will refer to E(x, y) as a conformal density, to underline that 

it does not correspond to a physical density of the membrane. 

We observe that, if S(x, y) is a conformal density, its integral over the square is just the area of the region T>, 



For s — > 1 + we may obtain a general behavior of the spectral zeta function to order O [s — 1) as 

n\T,\n) 



E 

n=l 

oo 



+ 0(s-l) 

(nlSln) 



E 

n=l 



4L 2 e r , 



•O(a-l) 



where the second series converges for s = 1. The first series, on the other hand, reproduces the spectral zeta function 

'A 



of a square of side 2L, apart for a multiplicative factor (jj?)" ■ For this reason, we may obtain the general behavior 



of the spectral zeta function at s = 1 as: 



Z(s) 



where 



9 = 9n + log 



A An 
AL 1 + ^A 



E 

71=1 



A_ 

Air 



1 



s - 1 



.4 
4L 7 



, M V 7i (|) -7i(i) 27T 

log U^) + * T 



- + 5e 



(nlSln) 



.4 

4L 7 



Itzykson, Moussa and Luck have found this behavior long time ago ( eq.(6) of ref. Q), obtaining an expression for 
g in terms of an integral. 
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3. Spectral zeta function of a deformed squared membrane 



As a specific example, we consider the domain obtained from the map 



/(*) = 



z + az 



1 + 



8a 2 L 2 



for \a\ < 1/2 (we have studied this problem for \a\ < 1 in a previous paper, ref.[30[, calculating the lowest energies 
of this domain in perturbation theory and comparing them with the numerical values obtained using a collocation 
approach). The domain corresponding to a = 1/2 is shown in Figf4] 

The map is defined to allow A = An ; the conformal density in this case is 



S(x,y) 



3 (4c*V + (2ax + l) 2 ) 



8a 2 L 2 + 3 



and its matrix elements in the basis of the square are 

1 



(n x ,n y \T,\n' x ,n) 



6Aa 2 (nl+n 2 y ) 
7r 2 nln 2 {2Aa 2 +3) 
48Aa 2 n v n' y (^(-l) n y + n y +l) 
7r 2 (2A a 2 +3)(n v -n y ) 2 (n y +n y ) 2 
48v^4anxn^(v / ^a((-l)™ x + "^+l)+(- 







(2Aa 2 +3) (n x -n'J 2 (n x +< ) 2 



n x 






n y 




n y 


n x 




K 


n y 






n x 




K 


fly 




n' y 






K 


n y 




n' y 



Using the diagonal matrix elements above, one finds in this case 

A\ 71 (I)- 71 (7) 2n 



log 



4tt 2 



+ 7 



2nAa 2 



3 (2^q 2 + 3) 



In Table U we report the values of the spectral zeta function at s — 2 for the deformed square for different values of 
a: the exact value, Z(2), is compared with the value obtained considering only the diagonal contributions, Z{2)^ dlaa \ 
and with the values obtained from the first 3000 numerical eigenvalues obtained with a collocation approach with a 
grid with 9801 points (adding the contribution of the higher excited states via the Weyl's law): 

2 

n=3001 
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where A — 4 is the area of the domain and C its perimeter; in the last column we report the value obtained only 
using Weyl's law: 

oo 

Z^Weyl) = 

n=l 

Notice that Z(2)^ dla9 " > accounts almost completely for the value of Z(2), up to a w 1/10. 



a 


Z(2) 




Z{^2^ { nurn ) 


Z(2) {Weyl) 


i 

100 


0.06970508 


0.06968939 


0.06970508 


0.04950760 


i 

50 


0.06969987 


0.06963720 


0.06969987 


0.04950758 


1 

25 


0.06967869 


0.06942953 


0.06967869 


0.04950735 


1 

10 


0.06951485 


0.06802167 


0.06951486 


0.04949801 


1 

4 


0.06805735 


0.06073539 


0.06805740 


0.04918833 


1 
2 


0.06143122 


0.04641541 


0.06143131 


0.04662053 



TABLE I. Spectral zeta function at s — 2 for the deformed square of area A — 4 for different values of a. 




4- Spectral zeta function of the annulus 

The function 



m = e* 



-L 



maps the rectangle [— L,L] x [— ir, ir] onto a circular annulus of external radius R = 1 and internal radius r = e 2L . 
The conformal density in this case is 



,y) 



and depends only on x. 

Since we are interested in calculating the Casimir energy of this configuration we need to evaluate the contributions 
of both the transverse electric (TE) and transverse magnetic (TM) modes. The first ones correspond to imposing 
Dirichlet boundary conditions at the border of the annulus, while the second ones correspond to imposing Neumann 
boundary conditions at the border. 

In the case of TE modes we use the basis on the rectangle is (see ref . [40|] ) 



Xn y (y) , u = 1 
<Pn y (y) , U = 2 



^n % ,n y ,u{x, y) = 1pn x ( x ) x 

where 

[ / IX 7T \ 



and 




* 9 (28) 



1,2,.. 



which corresponds to Dirichlet boundary conditions at x = ±L and periodic boundary conditions at y — ±vr. Each 
state is determined by three quantum numbers, n x , n y and u. Alternatively we may also impose Dirichlet boundary 
condition on both directions: this case corresponds to a circular annulus with a cut. 
The eigenvalues of the negative laplacian on this domain are 



(DP) _ , 2 n 
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and 

(DD) n 2 x TT 2 n 2 y 

where the superscripts (DP) and (DD) refer to Dirichlet-Periodic and Dirichlet-Dirichlet boundary conditions respec- 
tively. We appreciate that the states DP with n y > 1 are doubly degenerate, whereas states with n y = are non 
degenerate. 

The matrix elements of £ are the same in both basis since the density does not depend on y: 

(n x ,n y) u\E\n x ,n y ,u')( DP > = (n x , n v \E\n' x , n y ) (DD) 

7r 2 n 2 (r 2 -l) _ , _ , _ , 

2(7r 2 n 2 log(r)+log ;i (r)) ' n x — n x ,n y — U yl U - U 
S^n^logCr) ((-!)"* + < -r 2 ) _ ; _ f ' 

"(7r 2 (n x -n^) 2 +4 log 2 (r))(7r 2 (n, +<) 2 +4 log 2 (r)) ' ^ ^ ^ ' n ^ _ ' M _ U 

Notice that the (DD) matrix elements only hold for n y > 1. 

We will first concentrate on the circular annulus (DP): the value of its spectral zeta function s — 2 is 



Z^ DP ){2) — y-v (n x , n y , u\T,\n' x , n y , u)' 



(DP) (DP) 
n^,n y ,u n' x tn xi n yi ut n' :1 .,ny,u 

In the case of the circular annulus the exact eigenfunctions of the negative laplacian are known and read 



cos tlO , u = 1 
sinn# , u = 2 



where r < p < 1; N mns is a normalization constant and J and V are the Bessel functions of first and second kind. 
The corresponding eigenvalues are obtained from the solutions of the equation 

Y m (k) J m (kr) - J m (k)Y(kr) = , (29) 

as 

P — b 2 
J^mn — n-mn • 

In ref.[40| we have observed that the energies of thin annuli are approximated extremely well by the formula 

ADP) 



(n x ,n y ,u\T,\n x ,n y ,u} 
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This formula not only contains the correct asymptotic behavior (Weyl's law) but it also describes quite accurately 
the low energy behavior of the eigenvalues (see for instance Fig. 2 of ref.[4(|). In this case one expects that most of 
the contribution to Z(2) must be carried by the diagonal terms: 

z K 9J (2)= 2^ ; (dp) r 2 • ( 3 °) 

n x ,n B ,u ,n v ,u) 

In TablelHwe compare the exact value Z i - DP ^{2), with Z(*° s >(2) and with 




where A = 7r(l — r 2 ) is the area of the circular annulus and L = 2-7r(l + r) its perimeter. En are the first 10000 
numerical eigenvalues obtained solving the trascendental equation (|29p . each calculated with a precision of 20 digits. 



r 


Z(2) 




Z{^2^ ( nurn ) 


Z(2) (Weyl) 


1/10 
1/2 
9/10 


0.0257710759 
0.0057419570 
0.0000578599 


0.0169570674 
0.0054705758 
0.0000577934 


0.0257710743 
0.0057419569 
0.0000578601 


0.030450016 
0.011541075 
0.000251957 



TABLE II. Spectral zeta function at s — 2 for a circular annulus of varying internal radius r and fixed external radius R = 1. 



Notice that for r — 9/10 the off-diagonal contributions account only for about 1/1000 of the value of Z(2), confirming 
the observation made in ref. 40]; the value obtained using Weyl's law on the other hand largely overestimates Z(2), 
since this approximation fails to describe the low energy part of the spectrum. 

It is interesting to obtain the quantity 



AZ(1)= Urn Z^ DP \s) - Z^ DD \s) 

S— >1+ - 



which is finite. 
We have 



oo oo 



A<7m o,i|£K,o,i} 2 , . . , s 



n w = l \ 2L I n w = l n y =l 

(1 - e~ 4i ) (2Lcoth(2L) - 1) 



L 4L 2 



V 4L 2 4-1 



E 



i x sinh(2L)csch (^A 



1 



81og(r 
1 
2 



16L 

[(l-r 2 ) + (l + r 2 )log(r)] 



4L 2 + 7r 2 n 2 



5(1- OS 

3=0 



B e , l - 



i log(r) 



) +B\ e EOT , 



cy+H* 3 *log(r) 



1,0 



where B(z,a,b) is the incomplete beta function. This expression is exact and holds for all values r, < r < 1: in 
particular, around r = 1 we have 

AZ(l) « — (l-r) 2 . 

We will now discuss the analytic continuation of the spectral zeta function of the annulus, in the limit r — > 1~ 
(L — > + ). In this limit, it is possible to approximate the conformal density as 



S(ar,i/) w(l-2i)(l + 2a?) + 
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and 

SE(x, y) w -2L + 2x(l - 2L) + . . 

Working to this order we have that 

(n x ,n y ,u\8'E\n x ,ny,u) = -2L 
and the spectral zeta function of thin annulus may be approximated as 



1 - 2sL ^ 1 - 2sL ^ ~ l _ 2sL 



z^ dp \s)^ y ; s = y ; +2V y 



^ ( (DP) \ ^ ( (DP)V ^— ' ^— ' / (DP) V 



An analogous expression can be obtained for the case of a thin annulus with a trasversal cut, which corresponds to 
choosing Dirichlet boundary conditions along the two directions. In this case we have 



(DD) _ "a 71 " 2 , n v 



where the superscript (DD) indicates Dirichlet-Dirichlet boundary condintions, and n Xl n y = 1,2,.... The spectral 
zeta function in this case reads 



1 - 2sL A A 1 - 2sL 



r (dd) y ^ ^ / (dd) \ s ' 

The analytic continuation of these expressions is discussed explicitly by Kirsten in ref . [4l[ : following [Zlj we define 



£ 2 = 1£ 3 = 1 \ \ ■* / \ J / 



The analytic continuation of this function is given in eq.(4.31) of |41| : 

Cc (L 2 ,L 3 , S ) = --^-j C(2«)+ 2V _ r(s) ^-j C(2— 1) 



Therefore 



Z^ (s) « (1 - 2sL) [4 s tt- 2s j L 2s C(2s) + 2 Cc (tt, 2L, s)] 
z( DD )( S )«(l-2.sL)Cc(2^,2L, S ) . 
Around s = —1/2, in the limit r — > 1 — , we have 

Z^ p \-l/2) « z(° D )(-l/2) « + • • ■ , (31) 

which provides the transverse electric contribution to the Casimir energy of an arbitrarily thin annulus: 

„(TE) _1 7 {DP)( mo\~ C(3) r „„s 

It is easy to calculate the transverse magnetic contribution to the Casimir energy keeping in mind that e„i = e^°^ 
and f (NN) - f (DD) 

4™) = i^)(-l/2 ) «-^ F + ... (33) 
Thus the electromagnetic Casimir energy of a thin annulus is 
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Spectral zeta function of concentric cylinders 



Gosdzinsky and Romeo [42| have related the spectral zeta function of a circle to the spectral zeta function of a 
cylinder of circular section; Nesterenko and Pirozhenko [43| have also used this formula to obtain the Casimir energy 
of a circle using the results for a cylinder. 

Using our conventions, the relation between the two spectral zeta functions, eq.(3) of [HJ and eq.(6) of [43^ ) can be 
cast in the form 



Z cy i{s) 



1 

2^ 



B 







7 ( l \ 






^circle 1 « ^ J 



(35) 



Although the authors of refs. [42|, |43[ have applied this equation to the cylinder of circular section, the formula is 
general and it applies to cylinders of arbitrary section: therefore we can obtain the spectral zeta functions for these 
cases, working in two dimensions with a domain corresponding to the transversal section of the cylinder and then 
using eq. ([33|) . 

We write eq.(f3"51 in the general form as 



Zcyi(s) 



1 

2^ 



B 



1 



-, s - 



(36) 



We can apply eq. (1361) to calculate the Casimir energy of an infinite cylinder with annular section, in the limit where 
the two radii are very close (r — > R~ , R = 1). The spectral zeta function of the two dimensional domain has been 
already calculated to leading order in the previous example, and thus the calculation only requires the straightforward 
application of eq. (f31)| : 



ADP) 



,(DD) 



cyl I 2 I c V l 

The total electromagnetic Casimir energy will therefore be 



1 



360 (1 



ADP) 
J cyl 



7 {NP) 



360 (1 



(37) 



(38) 



The Casimir energy for concentric cylinders has been calculated by Mazzitelli and collaborators in 44]: our result 
reproduces eq.(61) of that papei0. This provides a useful check of the correctness of our approach. 

The same approach can also be applied to the calculation of the Casimir energy of cylinders and annular cylinders 
with slightly nonuniform density. 



V. CONCLUSIONS 



In this paper we have discussed a novel approach to the calculation of the spectral zeta function and of the heat kernel 
associated with the eigenvalues of a slighlty inhomogeneous system in a d-dimensional cube within perturbation theory. 
The method that we have devised relies on the standard Rayleigh-Schrodinger perturbation theory for hcrmitian 
operators and provides an explicit expression for the spectral zeta function associated to these operators that converges 
for s > d/2. To the best of our knowledge, this is the first calculation where the spectral zeta function of a system 
with arbitrary density perturbation has been obtained. 

The analytic continuation of this perturbative series allows one to obtain an approximation to the Casimir energy of 
the system under consideration. We have illustrated our method with some examples, evaluating the Casimir energy 
of different systems in one, two and three dimensions. In particular, for a string with piecewise constant densit y an d of 
two perfectly conduting, concentric cylinders of similar radius we have reproduced results already published [24 l44j. 
This provides a useful check of our approach. 

An interesting outcome of our analysis is the appearance of irremovable divergences already to first order pertur- 
bation theory in the Casimir energy of one dimensional inhomogeneous systems, when either Dirichlet-Dirichlet or 
Neumann-Neumann boundary conditions are applied at the ends of the string, unless the density £(&) of the string 
is such that £'(— L) = £'(£). The physical interpretation of these singularities is not clear. On the other hand, we 



5 Notice however that the calculation of [44| is not restricted to the particular limit that we are studying, r — > 1 
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also have found that the sum of Dirichlet-Dirichlet and Neumann-Neumann zeta functions is free of divergences, to 
first order. 

For positive integer values of s, with s > d/2, our formalism provides exact (non-perturbative) sum rules that can 
be useful to test numerical and analytic approximations to the eigenvalues of the problem. 

The calculations presented in this paper have been mostly performed to first order in perturbation theory: the 
analytic continuation of the higher order expressions, which involve multiple series, is a non trivial task for general 
densities (in the paper we have discussed a one dimensional problem where the second order calculation is easily 
obtained) Q We are currently studying a non-perturbative extension of the present approach that may allow to 
describe a larger class of problems. 
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